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ABSTRACT 

We present HST STIS observations of two extreme NLSls, IRAS 13224-3809 and 
IH 0707—495. The spectra are characterized by very blue continua, broad, strongly 
blueshifted high- ionization lines (including C IV and N V), and narrow, symmetric 
intermediate- (including C III], Si III], Al III) and low- ionization (e.g., Mg II) lines 
centered at their rest wavelengths. The emission-line profiles suggest that the high- 
ionization lines are produced in a wind, and the intermediate- and low-ionization lines 
are produced in low-velocity gas associated with the accretion disk or base of the wind. 
In this paper, we present the analysis of the spectra from these two objects; in a com- 
panion paper we present photoionization analysis and a toy dynamical model for the 
wind. The highly asymmetric profile of C IV suggests that it is dominated by emission 
from the wind, so we develop a template for the wind from the C IV line. We model 
the bright emission lines in the spectra using a combination of this template, and a 
narrow, symmetric line centered at the rest wavelength. We also analyzed a comparison 
sample of HST spectia from 14 additional NLSls, and construct a correlation matrix of 
emission line and continuum properties. A number of strong correlations were observed, 
including several involving the asymmetry of the C IV line. 

Subject headings: quasars: emission lines — quasars: individual (IRAS 13224—3809, 
IH 0707-495) 
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1. Introduction 

In 1992, it was demonstrated by Boroson & Green that the optical emission hne properties 
in the region of the spectrum around H/? are strongly correlated with one another. A principal 
components analysis allowed the largest differences among optical emission line properties to be 
gathered together in a construct commonly known as "Eigenvector 1". The strongest differences 
hinged on the strength of the Fc II and [O III] emission, and the width and asymmetry of H/3. 
This set of strong correlations is remarkable, as it involves correlations among the dynamics and 
gas properties between emission regions separated by vast distances. Furthermore, the same set 
of correlations are observed in samples selected in many different ways. This pervasiveness leads 
us to believe that we are observing the manifestation of a primary physical parameter. A favored 
explanation is that it is the accretion rate relative to the black hole mass onto the active nucleus 
(e.g.. Boiler, Brandt & Fink 1996; Pounds, Done & Osborne 1995; Comastri et al. 1998; Leighly 
1999ab; Mathur 2000; Mineshige et al. 2000; Kuraszkiewicz et al. 2000; Pounds et al. 2001; 
Puchnarewicz et al. 2001; Ballantyne, Iwasawa, &: Fabian 2001; Boroson 2002; Crenshaw, Kraemer 
& Gable 2003). 

Narrow- line Seyfert 1 galaxies (NLSls^) are identified by their optical emission line properties. 
They typically have narrow permitted optical lines (FWHM of H/3 < 2000 km s~^), weak forbidden 
lines ([O III]/H/? < 3; this distinguishes them from Seyfert 2 galaxies), and frequently they show 
strong Fe II emission (Ostcrbrock & Pogge 1985; Goodrich 1989). These arc the same properties 
that define the Boroson & Green (1992) Eigenvector 1, and indeed, NLSls fall at one end of 
Eigenvector 1. Thus, the study of NLSls offers an attractive research opportunity: if we can 
understand the origin of the emission-line properties of NLSls, then we may be in a position to 
understand AGN emission in general. Furthermore, because as least some of the lines in NLSls 
are narrow, identification of the frequently strongly-blended lines is less ambiguous than it is in 
broad-line objects. 

NLSl research dramatically increased after the publication of a seminal paper in 1996 by Boiler, 

Brandt & Fink describing their discovery that NLSls have systematically steeper ROSAT spectia 
than AGN with broad optical lines. Subsequently it was shown from ASCA observations that the 
hard X-ray spectrum is also steeper (Brandt, Mathur & Elvis 1997; Leighly 1999b), and that the 
X-ray emission shows higher fractional amplitude of variability (Leighly 1999a). These results were 
generally interpreted as evidence that the black hole is systematically smaller in NLSls, implying a 
higher accretion rate relative to the Eddington value compared with Seyfert 1 galaxies with broad 
optical lines. 

Once the high accretion rate paradigm was established, it could be used to understand why 
NLSls have narrow optical lines. Boiler, Brandt, & Fink (1996) mention that if the black hole mass 



^In this paper, the term NLSls is used to refer to aU objects fulfiUing the optical criteria regardless of their 
luminosity. Likewise, the term broad-line Seyfert 1 is used to refer to objects of both Seyfert and quasar luminosities. 
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in NLSls is smaller, and the emission lines are produced at the same radius as in broad-line objects 
(due to the same radiating flux) the dynamical velocities will be smaller, and narrower lines will be 
produced. This theme was expanded upon by Wandcl &: Boiler (1998), who consider the fact that 
the steeper soft X-ray spectrum of NLSls, when extrapolated into the UV, implies a larger ionizing 
flux relative to their luminosity. This causes the emission region to be larger where the velocities 
are smaller. 

However, the optical emission lines are only one side of the coin; UV emission lines, generally 
characterized by a higher degree of ionization, may also provide valuable information. This is 
because the optical and UV emission lines are generally not produced in the same gas. The best 
evidence for this comes from so-called reverberation mapping experiments (e.g., Peterson 1993). 
Also, the profiles of the emission lines suggest this as well, as the higher ionization lines are observed 
to be generally broader than the lower ionization lines (e.g., Baldwin 1997). Therefore, we may 
learn more about NLSls by studying their UV emission lines. 

Several groups have investigated the UV emission line properties of NLSls. The first paper 
investigating a sample of lUE spectra from NLSls was published by Rodriguez-Pascual, Mas-Hesse 
& Santos-Lleo (1997). They reported that the strongest UV emission lines have broad wings, unlike 
the optical emission lines, and they speculate that the UV line-emitting gas may be optically thin 
in the continuum. 

Kuraszkiewicz et al. (2000) present a detailed study of lUE and HST spectra from NLSls. 
They found that, compared with broad-linc quasars, the equivalent widths of C IV and Mg II are 
low. They found that, in the cases where the lines could be dcblcndcd, the Si ni]/C III] ratio was 
high. In addition, the Si IV-l-0 IV] feature near 1400 A was found to be strong compared with 
C IV, and the Al III equivalent width was high. They also found that although the UV emission 
lines in NLSls are broader than the Balmer lines, they still tend to be narrower than those lines in 
broad-line quasars. They then model many of the emission line ratios using a 1-zone Cloudy model 
(Ferland 2001). They infer that the line ratios imply that the ionization parameter is lower and the 
density is higher in NLSls than in broad-line Seyfert 1 galaxies. We comment on their inferences 
in the companion paper to this one (Leighly 2004a; hereafter referred to as Paper II). 

Wills et al. (1999), studying the HST spectra from a small sample of PG quasars, found that 
Eigenvector 1 is also associated with a high Si III] to C III] ratio, stronger UV low-ionization lines, 
weaker C IV, but stronger N V. They infer from these spectra that NLSls are typified by high 
densities, low ionization parameters, and nitrogen enhancement, perhaps from nuclear starbursts. 

All the work done so far on the UV spectra of NLSls has been generally confined to studies 
of the properties of the integrated lines, to a greater or lesser extent, depending on the analysis. 
Specifically, the conditions of a single-zone BLR are modelled or discussed more or less quanti- 
tatively, and evidence for a radially-extended BLR is inferred when particular lines arc found to 
not be consistent with the one- zone model (e.g., Kuraszkiewicz et al. (2000) infer that, since the 
observed Mg II is underpredicted by the one-zone model, it must be produced in a different region 



-4- 



than the other UV emission hnes). However, we know that generally the emission line profiles are 

not all the same, and higher ionization lines tend to have broader profiles. This implies that the gas 
properties are different among regions with larger or smaller line-of-sight velocities, and therefore 
a study of the average properties could be misleading. 

In this paper, we present the observations and analysis of HST observations of two NLSls 
that have truly remarkable UV spectra. In Section 2, we introduce the targets, describe our 
observations and discuss the analysis of the emission lines. In Section 3, we compare the results 
from these two NLSls with spectra from 14 NLSls drawn from the HST archive, and perform a 
correlation analysis and a principal components analysis on the results. In a companion paper II, 
photoionization analysis, and a toy dynamical model are presented, and the results are discussed. 
Partial results from this study have been published previously in Lcighly (2000) and Leighly (2001). 
We use Hq = 50kms~^ Mpc~^ and qo = 0.5 unless otherwise noted. 

2. Observations and Analysis 
2.1. The Targets 

We observed two NLSls, IRAS 13224-3809 (z=0.066, V=15.9) and IH 0707-495 (z=0.040, 
V=15.7), using HST STIS during 1999. Groundbased optical spectra were also obtained. The 
observing log is given in Table 1. 

These two NLSls are quite remarkable objects and can be considered to be extreme for their 
class in terms of their X-ray properties. They showed among the largest fractional amplitudes of 
variability in a sample of NLSls observed by ASCA (Leighly 1999a). i? 05^ T monitoring observa- 
tions of IRAS 13224-3809 found strong soft X-ray flaring by up to a factor of 60 in as short time 
interval as two days (Boiler et al. 1997). Large amplitude X-ray variability was also reported dur- 
ing a long ASCA observation (Dewangen et al. 2002). Their X-ray spectra are also unusual. They 
show extremely strong X-ray soft excess components and lie on the extreme of the variability-soft 
excess correlation (Leighly 1999b; Leighly 2004b). XMM-Newton observations of IH 0707—495 
and IRAS 13224—3809 reveal a peculiar deep absorption feature near 7 keV that is quite difficult 
to explain (Boiler et al. 2002; Boiler et al. 2003). The XMM-Newton observation of IH 0707-495 
showed primarily color-independent X-ray variability (Boiler et al. 2002). In contrast, a Chandra 
observation of IH 0707—495 found strong spectral variability in that the soft X-rays varied much 
less than the hard X-rays (Leighly et al. 2002; Leighly et al. in prep.). 

The motivation for the HST observations was the discovery of a peculiar feature near 1 keV 
observed in the ASCA spectra (Leighly et al. 1997). This feature appeared to be an absorption edge, 
but the energy was clearly too high to be ascribed to highly-ionized oxygen, which is generally the 
most important source of X-ray opacity. Leighly et al. (1997) speculated that oxygen was indeed 
responsible for the feature, and, if so, the absorbing gas was outflowing at rather large velocities. 
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This interpretation was strengthened by the similarity between luminous NLSls and low-ionization 
BALQSOs (e.g., Boroson & Meyers 1992 for I Zw 1). Reasoning that X-ray absorption may be 
accompanied by UV absorption, we proposed to make use of the much superior spectral resolution 
(relative to ASCA CCDs) provided by HST STIS, and look for evidence for outflowing gas in the 
UV. We note here that no evidence for absorption intrinsic to the quasars was found in the spectra. 
This result leaves the issue of the X-ray features unresolved, since it is certainly plausible that the 
X-ray absorption lines could be produced in gas too highly ionized to present opacity in the UV. 
We also note that subsequent high resolution X-ray observations using Chandra reveal absorption 
lines (Leighly et al. 2002; Leighly et al. in prep.) 



2.2. The HST observations and Continua 

Our observations were made during the first cycle that the STIS spectrometer was available. 
We made the HST observations using a "pattern" in order to decrease possible systematic error 
due to fixed pattern noise. Three or four steps along the slit were used for each spectroscopic 
element. The relative wavelength calibration of each spectrum was checked prior to summing the 
spectra by cross correlating with a preliminary average spectrum. Many Galactic absorption lines 
were observed in the spectra. These were used to identify and correct zero-point wavelength offsets 
in the spectra. Then, the redshift was obtained from the optical and UV spectra. In the optical 
spectrum of IRAS 13224—3809, H/? is extended in the spatial direction. The extended emission 
probably originates in a starburst in the host galaxy, and therefore provides a good estimate of 
the systemic redshift. The extended emission contributes to the narrowest part of the H/3 profile; 
therefore, we estimated the redshift from centroid of the top of the II/3 line. Interestingly, O III 
has a blue offset of 370 kms^^ (Grupe &: Leighly 2002), similar to I Zw 1 (Phillips et al. 1976). 
The redshift of IH 0707—495 was estimated also using the centroid of the top of the H/3 line. 
These redshifts were consistent with that inferred from Mg II, if that line is centered at the rest 
wavelength, and the centroid of the narrow part of Lya. The measured redshifts were 0.066 and 
0.040 for IRAS 13224-3809 and IH 0707-495, respectively. 

The UV spectra are displayed in Fig. 1 along with nonsimultaneous optical spectra obtained 
at CTIO (details in Table 1) and a composite quasar spectrum from the LBQS (Francis et al. 1991) 
for comparison. This figure shows that IRAS 13224—3809 and IH 0707—495 have continua as blue 
as the average quasar. This result is in agreement with the discovery by Grupe et al. (1998) that 
soft X-ray selected AGN from the ROSAT All Sky Survey, a sample which is comprised of 1/2 
NLSls, generally have quite blue optical-UV continuum spectra; some are as blue as predicted 
by the standard (thin) accretion disk model {F^, oc i/ ' , Fx (X X-^/^). For comparison, a line 
with a slope of —7/3 is included on Fig. 1 (which is plotted in terms of Fx), as well as one also 
with slope —0.32, which is the median slope in the LBQS (Francis et al. 1991). Very blue optical 
spectra seem to be common among luminous NLSls (see also RGB J0044-|~193: Siebert et al. 1999; 
RX J0134.2-4258: Grupe et al. 2000; RX J2217.9-5941: Grupe, Thomas & Leighly 2001; PHL 
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1811: Leighly et al. 2001). 

The optical spectrum from IH 0707—495 appears to be slightly offset in normalization from 
the UV spectra. This offset could be due to aperture effects, or from variability, as there was a one- 
year interval between the optical and UV spectra. However, the continuum slopes in the UV and 
optical spectra appear to be consistent. In contrast, the optical spectrum from IRAS 13224—3809 
appears to be much flatter than the UV continuum spectrum. This object is an extremely luminous 
far-infrared source (Norris et al. 1990) , and we suspect that the flattening is due to contamination 
of the central engine emission by the star burst. We suspect that similar contamination was present 
in other spectra from this object, leading to reports that the optical-UV continuum spectrum is flat 
(e.g., Boiler et al. 1993). The UV spectra presented here show that the central engine continuum 
spectrum is not flat, but very clearly blue. 

2.3. The Emission Lines and Modeling 

The spectra are shown in more detail in Fig. 2. In this plot, we follow Laor et al. (1997) and 
mark the rest wavelengths of many of the strong lines frequently observed in AGN spectra. Our 
identification of the emission features are discussed below. 

These spectra are characterized by the following properties: 

• Strong Fe II emission. Emission lines from Fe III are also identified. The Fe II is strong 
enough that it forms a pseudocontinuum makes identification of the real continuum difficult, 
but it can also be identified as a forest of relatively narrow lines, especially near 2000 A. 

• Broad and blueshifted high-ionization lines. The high-ionization lines, including Lya, N V, 
the Si IV-0 IV] complex, and C IV are all relatively broad (FWHM are approximately 
5000 km s~^) and blueshifted, with little emission redward of the rest wavelength, at least in 
C IV and N V. 

• Relatively narrow intermediate- and low-ionization lines centered at the rest wavelength. 
C III], Si III] and Mg II are narrow (900-1900 km s"^), and their profiles are reflection- 
symmetric about about the rest wavelength. 

• Strong lower ionization lines, including Si II and C II. 

The spectra most strongly resemble that of the NLSl prototype I Zw 1, which also shows 
blueshifted high-ionization lines, narrow intermediate- and low-ionization lines and strong Fe II 
and Fe III (Laor et al. 1997). However, as discussed in Section 3, they are somewhat different than 
the other narrow-line Seyfert 1 galaxies and quasars. The presence of strong low-ionization lines in 
our spectra is similar to that seen in several PG quasars (Wills et al. 1999; see their Figure 1). 
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The profound differences in the profiles of the high- and low-ionization hnes are highhghted 
in Fig. 3, which compares the profile of the prototypical high-ionization line C IV with the profile 
of a low-ionization line Mg II, for our NLSls and for a composite quasar spectrum^ (Francis et al. 
1991). It has been previously noted that high-ionization lines are frequently broader and strongly 
blueshifted compared with the low-ionization lines (e.g., Tytler & Fan 1992; Marziani et al. 1996). 
The velocity offsets in C IV described by Marziani et al. (1996) range from zero to ~ 3000; the 
centroids of our profiles are offset by ^ 2500 km s~^, as large as the largest values values found 
by Marziani ct al. These profiles are most simply explained if the high-ionization line-emitting 
gas is accelerated toward us in a wind, while the low-ionization line-emitting gas is located on the 
surface of the accretion disk, or in the low-velocity base of the wind. The wind may be accelerated 
by radiative-line driving from the strong UV continuum. This process is commonly inferred to be 
occurring in hot stars and cataclysmic variables, but it has also been applied to AGN (e.g., Murray 
et al. 1995; Murray & Chiang 1998; Proga, Stone, & Kallman 2000), and we discuss this further 
in Paper II. The accretion disk is assumed optically thick, so that we do not see the receding wind 
on the other side. Hence, the emission lines arc predominantly blueshifted. Collin- Souffr in ct al. 
(1988) were the first to suggest a two-component broad-line region in which the high-ionization 
lines are produced in a wind, and the low-ionization lines are associated with the outer parts of the 
accretion disk. The spectra presented here arguably provide the strongest support for this scenario 
to date. 

Because the high-ionization lines are nearly disjoint in velocity space compared with the low- 
and intermediate- ionization lines, we model them separately. We then use the results to constrain 
photoionization models for the conditions of the disk and wind separately in Paper 11^. Such 
an undertaking is fraught with peril. In the emission-line profile, we observe only the velocity 
component parallel to our line of sight. Especially problematic is the gas at zero velocity, as 
the emitting matter may have huge transverse velocity that we cannot detect. However, such an 
approach is not without precedent: Baldwin et al. (1996) analyzed the spectra from a luminous 
quasar Q0207— 398 in this way. 

To measure the line properties, we adopted a modified version of the procedure used by Baldwin 
et al. (1996). Wc outline the procedure here, and describe the application to each region of the 
spectra below; the results are given in Fig. 4 and Table 2. First, we make several simplifying 
assumptions. Fig. 3 suggests that, to zcroth order, the wind and disk lines are kinematically 
disjoint. We further assume, to begin with, that all high-ionization lines have the same velocity 



■^Note that composite spectra constructed on the basis of redshifts measured from a single hne may have broader- 
than-average and more-symmetric-than-average hnes. 

'Wc do not know with certainty the geometrical and physical origin of the emission lines in the objects we arc 
discussing here. However, for simplicity, we refer to the highly blueshifted high-ionization lines £is originating in the 
"wind" , and the narrow, symmetric low-ionization lines as originating in the "disk" . These distinctions are somewhat 
similar to the HIL and LIL regions previously proposed (Collin-SoufFrin ct al. 1988), except that C III] and other 
similar intermediate-ionization lines in our spectra appear to be produced in the disk. 
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profile; we comment on the validity of that assumption below. Then, since C IV is isolated compared 
with the other high-ionization hncs, wc use it to create a wind-profile template that is used to model 
the other high-ionization lines. The procedure used to develop the template is described in the next 
section. 

The intermediate- and low-ionization lines appear to be predominately symmetric about their 
rest wavelengths; therefore, we model them using symmetric Lorenzian or Gaussian profiles. Pre- 
liminary examination of spectra showed that we could not assume the same velocity width for all 
the intermediate- and low-ionization lines, as Mg II could be seen to be clearly significantly nar- 
rower than C III]. Therefore, in fitting the intermediate- and low-ionization lines, we allowed the 
width to remain an adjustable parameter so long as it could be constrained usefully in the fitting. 
Generally, the widths could be usefully constrained unless the lines are very weak (e.g., Section 
2.3.7), or unless the feature is comprised of several heavily-blended components, in which case we 
fixed the velocity width for particular individual lines to the values measured from other lines. 

The next sections explain the process that we used to model the spectra. Throughout, we 
use a locally-defined continuum. This is necessary because the strong Fe II and Fe III multiplets 
that may form a pseudo-continuum that masks the real continuum. We note that placement of the 
continuum is a potentially significant source of uncertainty in the measurement of the emission-line 
fiux. 

2.3.1. The CIV region 

First, the continuum was locally identified around the C IV emission line, then fit with a linear 
model and subtracted. Next, we removed the portions of the profile contaminated by absorption 
lines originating in our Galaxy (Blades et al. 1988). The Galactic column is significant in the 
direction of our quasars, and there are a number of fairly strong lines apparent in the spectrum. 
Careful examination of the C IV profile reveals faint absorption lines that were detected by Blades 
et al. but were not identified. The regions of the profile contaminated by absorption lines were 
removed and the profile was fitted with a high-order spline model. 

The C IV line is a doublet, but the final template should consist of the contribution of only one 
component. To derive the profile of one component, we used a procedure similar to that used by 
Baldwin et al. 1996. We assumed that each of the two components contributed equally to the line 
(that is, the line is optically thick). The approximate contribution of one component of the doublet 
was then obtained from the spline model by a sort of bootstrap method. The doublet components 
are separated by 4-5 wavelength bins. Therefore, the 4 or 5 points at longer wavelength side of the 
spline model consists solely of the lower energy component, so these can be used to determine the 
amplitude of the 4 or 5 longest wavelength points of the template. Points shortward of this consist 
of the sum of both lower-energy and the higher-energy components. At the next point, where i = 6, 
for example, the amplitude of the higher-energy component at i is already known, because it must 



-9- 



have the amplitude of the lower energy component at i =1-2. Then, since the sum of the lower and 
higher-energy components at any point i must equal the spline model at that point, the amplitude 
of lowcr-cncrgy profile at point i can be solved for. Thus, the amplitudes of the remaining points 
could be determined iteratively. The C IV emission lines, the contribution of each multiplet and 
the resulting profiles are shown in Fig. 4. Interestingly, the shape of the template is slightly but 
distinguishably different between the two objects. 

2.3.2. The Lya and N V region 

The template developed using the C IV profile, shown in Fig. 4, was next fit to the Lya/N V 
region in both spectra. A second-order polynomial was used to model the local continuum, which 
was then subtracted. The fitting was done by first transforming the template spectrum to the rest 
energy of the emission line in question and then adjusting the normalization until some part of the 
profile did not obviously exceed the height of the spectrum; i.e., if the model were subtracted from 
the data, no obvious negative residuals would be present. A narrow component of Lya was also 
clearly visible. This is inferred to originate from the disk, and it is modelled as a Gaussian. Both 
components of the N V doublet were modeled and the flux ratio between the doublet components 
was again assumed to be 1 to 1. The results are shown in Fig. 4. 

This procedure resulted in a good fit for IRAS 13224—3809. Excess emission shortward of Lya 
may come from C III* A1176 line, a line that has been identified in the spectrum of the prototype 
NLSl I Zw 1 (Laor et al. 1997). Si II is strong elsewhere in the spectrum and may contribute near 
1190 Angstroms as well. Note that the deficit near 1222 A in the rest frame is caused by Galactic 
absorption (Fig. 2). 

The fit results for IH 0707—495 are less satisfactory. Excess emission is present both redward 
and blueward of the N V line. The blueward emission could also be present in IRAS 13224—3809 
but it may be hidden by the Galactic absorption line near 1222 A. This excess emission could be 
an indication that our assumption that all of the resonance fines have the same profile is to some 
degree incorrect; however, we cannot clearly see the blue side of that line because it is blended with 
the broad Lya line, so we cannot clearly evaluate whether this assumption is good or not. We plan 
to address this question in a future paper involving the O VI line observed in the FUSE spectrum 
of IH 0707-495. 

In both cases, since the wind profile fits N V well, we conclude there is no disk contribution to 
this line. Both a disk and a wind component are necessary, on the other hand, to model Lya. 
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2.3.3. The Si IV /O IV] region 

The feature near 1400A is clearly broad and asymmetric and at least part of it must be 
emitted in the wind. This is a complicated feature in AGN in general, because it may be comprised 
of emission from both O IV] and Si IV. Si IV is usually a lower-ionization line than the other lines 
produced in the wind. The disk contributes low-ionization lines (e.g., Mg II, as in Fig. 3); therefore 
it is possible that there is a contribution from Si IV from the disk in the 1400 A feature. To try 
to account for this possibility in a general way, our model for this region consists of components of 
wind emission from each of Si IV and O IV] that are modelled using the template developed from 
C IV as discussed in Section 2.3.1. We also include components of disk emission from each of Si IV 
and O IV] that are modelled as Gaussians, as discussed above. There are thus 4 components in the 
model. 

We assumed that the emission ratios of the two components of the Si IV doublet have equal flux 
(optically thick). For O IV] we assumed optically thin ratios. These depend partially on density 
and were found to be 0.011:0.225:0.458:0.082:0.225 for the A1397.2, A1399.8, A1401.2, A1404.8 & 
A1407.4 emission lines respectively (Flower k, Nussbaumer (1975) and Nussbaumer & Storey (1982) 
assuming a density of 10^ cm~^. 

We first determined the upper limits of the four components listed above by overlaying each 
one on the line profile separately and adjusting the amplitude until some part of the component 
was as large as the line; i.e., no large negative fit residuals were permitted. Those upper limits are 
listed in Table 2. 

Next, we investigated whether or not the entire feature could have been emitted in the wind; 
that is, is the disk component necessary? We overlayed the sum of the O IV] and Si IV wind 
components on the profile, adjusting the amplitudes as described above. While the shortest wave- 
length, blueshifted part of the feature could be modeled by the sum of the wind components from 
O IV] and Si IV, the longest wavelength part, near the rest energies of O IV] and Si IV, could not. 
We inferred that the excess emission near the rest energies of these features is emission from the 
disk. To model the disk component, emission line blends for Si IV and O IV] were constructed 
using Gaussians, assuming again that Si IV is optically thick and O IV] is optically thin, as above. 
The widths of these Gaussians could not be independently determined in the fit, so we assumed 
a FWHM of 7 Angstroms, corresponding to 1500 km s~^. This is close to the value found for the 
1900 A lines (Section 2.3.4). 

The 1400 A feature is very complicated. We show a model in Fig. 4 that appears to fit the 
profile, and the parameters for this model are listed in the Table 2; however, we do not claim that 
this model is unique. All we can really conclude is that emission from both the disk and the wind 
is necessary to explain the profile; the wind portion is needed to explain the shorter wavelength, 
blueshifted part, and the disk portion is needed to explain the rest-wavelength part. But we really 
cannot tell how much of each component originates in each of Si IV and O IV] . The photoionization 
modelling presented in Paper II is helpful for understanding this feature. 
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Adjacent to the 1400A feature, there is excess emission near 1420 A that could be another 
emission line. A possible candidate could be S IV] AA 1416.93, 1423.86 (Harper et al. 1999). 

2.3.4. The Al III/C III]/Si III] region 

In the wavelength range between 1800 and 2000 A, we expect to find emission from C III], 
Si III] and Al III, as well as possibly Fe III and Fe II. Identification of the continuum in this 
region is very difficult because of the strong Fe II and Fe III lines. In IRAS 13224—3809 we clearly 
identify fairly strong contributions from C III] , Si III] and Al III. These lines appear to be narrow, 
symmetric, strongly peaked, and centered at the emission line wavelength. Thus, we infer they are 
from the disk and that no measurable wind component is present in the 1900A feature. We use 
the spectral fitting package LINER (Pogge &: Owen 1993) to model the lines. Wc first attempt to 
model the emission lines as Gaussians; however, they are too strongly peaked for this to provide 
a good fit. Therefore, we use Lorenzian profiles and obtain a better fit, although we note that 
Lorenzian profiles may overestimate the line flux due to the large contributions in the wings that 
are difficult to distinguish from the continuum. However, there remains excess emission around 
1925 A that most probably originates in the Fe III 34 multiplet. This multiplet consist of three 
lines at 1895.46, 1914.06, and 1926.30 A (e.g., Graham, Clowes & Campusano 1996). Wc include 
these three emission lines in the model, assuming that their fluxes are equal (e.g., Hartig & Baldwin 
1986; see also Section 3). In the final model, we hold the emission line wavelength fixed at the rest 
wavelength, the width of C III] , Si III] and Al III are constrained to be equal, and the widths of the 
Fe III 34 multiplet components are constrained to be equal. 

The spectrum from IH 0707—495 appears to be more difficult to model well. There seems to 
be more Fe III present in the spectrum of this object. A completely satisfactory fit could not be 
obtained. 

2.3.5. Hell 

This line was very difficult to identify and measure. Following the example of I Zw 1, and the 
other high ionization lines, we expected it to be strongly blueshiftcd and broad; i.e., from the wind. 
Indeed, we see no evidence for a narrow emission line at 1640 A. We defined a local continuum 
from apparently line-free region, and measured the line using the wind template developed from 
the C IV profile, as was done above. This measurement is not very satisfactory as there are other 
emission features of the same amplitude, probably originating in Fe II, near the feature that we 
identify as He II. 
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2.3.6. Mg II region 

The spectrum around Mg II is contaminated by broad humps of Fe II emission. This makes 
measurement of the flux and width of Mg II difficult. Following several other authors (e.g., Corbin 
Sz Boroson 1996; Dietrich et al. 2002), we develop an Fe II template from the medium-resolution 
HST spectrum from I Zw 1 following the procedure used by Vestergaard & Wilkes (2001). We 
adjust the amplitude of the Fc II template so that it fits the Fc II feature in the spectra. We then 
subtract the Fc II component, and fit the Mg II line using the LINER package. It turns out we 
cannot constrain the doublet ratio as was done for I Zw 1 (Laor et al. 1997), so we require the 
doublet component fluxes to be equal. We cannot distinguish between Gaussian and Lorenzian line 
profiles; however, we note that these lines are significantly narrower than C III] and Si III]. 

2.3.7. Other low-ionization lines 

The fluxes of several weak low-ionization lines were required for photoionization modeling of 
the disk component in Paper 2. All of the lines considered in this section appeared to be centered 
at the rest wavelength, and symmetric, so we infer that they are entirely produced in the disk. 
They are sufficiently weak that wc could only measure their flux with any confidence; their widths 
appear to be consistent with the other low-ionization lines. The C II lines at 1335A and 2327A were 
modeled using Gaussians, employing LINER to do the fitting. C II A1335A was slightly blended 
with an unidentified line at 1345A (also seen in I Zw 1; see Laor et al. 1997) so these two lines were 
modelled simultaneously. This procedure worked weh for IRAS 13224-3809. IH 0707-495 posed 
a more difficult case, as there were several strong absorption lines cutting through the C IIA1335A 
lines. The C II A2327A was more difficult to model in both cases. It appeared to be perched on 
a broad, blended Fe II feature, and blended with Fe III. We were able to measure it with some 
confidence for IRAS 13224-3809, but did not obtain a measurement for IH 0707-495. We also 
measured the flux and equivalent width of N III]A1750A. This was difficult, as it is barely detectable 
in IRAS 13224—3809, so this measurement should be considered to be nearly an upper limit. It was 
more clearly detected in IH 0707—495, and we measure the flux and equivalent width by defining 
the local continuum on both sides of the feature and integrating over the small excess at 1750A. 

2.3.8. H/3 region 

We also obtained groundbased optical spectra from these two objects. Near H/3, the spectra in 
Narrow-line Seyfert 1 galaxies is frequently cluttered with multiplets of Fe II that make it difficult 
to measure the properties of H/3 and [O III] . This problem is overcome by subtraction of an Fe II 
template developed from a high signal-to-noise spectrum from the prototype NLSl I Zw 1 (e.g., 
Boroson & Green 1992). The Fe II subtracted spectra are displayed in Fig. 5. The H/3 line was 
fitted with a Lorenzian profile to obtain the width. 
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2.3.9. Comparison with Composite Spectra 

In Table 2, along with the emission-line measurements from IRAS 13224+3809 and IH 0707—495, 
we list measurements of emission lines from three composite spectra: the LBQS composite (Francis 
et al. 1991); the radio-quiet composite from ^fST spectra (Zheng et al. 1997); a composite spectrum 
from the FIRST Bright Quasar Survey (Brotherton et al. 2001). In all cases, the equivalent widths 
of the prominent features in our spectra are the same or smaller than the equivalent widths in the 
composites, except for N V, which has larger equivalent widths in our spectra. Zheng et al. attempt 
to deconvolve N V from Lya; in that case, Lya in our spectra are about one third, and N V are 
nearly double the equivalent width of that in the Zheng composite. The whole 1400A feature and 
1900A feature have about the same equivalent width as the composite. Zheng et al. (1997) attempt 
a deconvolution of the 1900A feature; our spectra have the same Al III equivalent widths, larger 
by a factor of two Si III] equivalent widths, and much smaller (by a factor of 3) C III] equivalent 
widths than the Zheng composite. The largest differences arc in C IV and Mg II, which have much 
smaller (one fifth to one half) the equivalent width of the composite spectra. These comparisons 
underline the fact that the UV emission lines, excepting N V, are relatively weak in NLSls. 

3. Comparison with Other NLSls 

The UV and optical properties of IRAS 13224-3809 and IH 0707-495, including their lu- 
minosity, their emission line profiles, and their equivalent widths and ratios, are almost identical. 
How do their properties compare with those of other NLSls? This is an important question for two 
reasons. First, their X-ray properties, while again virtually identical to one another, are markedly 
different from those of other NLSls. In a study of the X-ray properties of NLSls using ASCA 
data, Leighly (1999b) found that there is a correlation between the prominence or strength of the 
soft excess^ and the amplitude of the variability, with IRAS 13224-3809 and IH 0707-495 dis- 
playing among the strongest soft excesses and highest-amplitude variability in the sample. Also, 
these objects arc two of the three that possess an unusTial spectral feature near 1 keV that has 
been interpreted as highly blucsliifted absorption feature (Leighly et al. 1997) or as strong Fe-L 
absorption (Nicastro, Fiore & Matt 1999). We would like to know if the differences between IRAS 
13224-3809 and IH 0707-495 and other NLSls that are found in the X-ray band are also seen in 
the UV band. The second reason we are motivated to compare these two objects with other NLSls 
involves our interpretation of our analysis of the UV spectra that is discussed in detail in Paper II. 

To address these questions, we retrieved the spectra from a heterogeneous group of 14 NLSls 
from the HST archive, and analyzed them. Analysis of properties of NLSl UV spectra has been 
previously presented by several groups. However, reanalysis is necessary because we wish to in- 



^Thc soft excess (e.g., Turner & Pounds 1989) is the generic term for the excess of emission in the soft (0.1-2.0 
keV) X-ray band over the power law obtained in the 2-10 keV band. 
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vestigate both the emission hne profiles and their equivalent widths and ratios. The most detailed 
investigation to date was presented by Kuraszkiewicz et al. (2000). Our analysis differs from theirs 
in that they confined their analysis to the emission- line fluxes and ratios, while wc are interested 
also in the emission-line profiles, in particular that of C IV, as well. A comparison of the blueshift 
of C IV with the width of H/3 has been published by Marziani et al. (1996). However, they do not 
consider the fluxes and ratios of other important lines in the spectrum. 

The following objects, listed in order of right ascension, were investigated: Mrk 335, WPVS 007, 
I Zw 1, Ton S180, RX J0134.2-4258, PG 1001+054, RE 1034+39, PG 1115+407, PG 1211+143, 
PG 1404+226, PG 1402+261, Mrk 478, Mrk 493, and Ark 564. These objects represent a range of 
redshifts from 0.024 to 0.235, and span two orders of magnitude in monochromatic luminosity at 
2500A; the luminosities of IRAS 13224-3809 and IH 0707-495 arc in the middle of the sample. 
Almost all of these spectra were obtained using the FOS detector and therefore we reprocessed all 
of them, using the Post Operational Archives (POA) calibration^ when appropriate. Ton S180 was 
observed with the STIS detector, and we retrieved those data from the archive using "on-the-fly" 
reprocessing. 

We checked the wavelength calibration of each spectrum using the Galactic absorption lines. 
The spectra were then shifted up to ~ 1.5A as necessary. Our first estimate of the redshift was 
obtained from the spectra directly, by fitting the Mg II lines. We assume that Mg II lines are 
centered at the systemic redshift. This is justified by the similarity between Mg II and H/3 emission 
in AGNs (McLure & Jarvis 2002), and the fact that all of these objects have narrow H/3 and Mg II. 
These redshifts were in good agreement with those listed in NED. The spectra were then shifted 
into the rest frame. The Galactic E(B-V) was obtained from NED, and the spectra were corrected 
using the COM reddening law (Cardelh, Clayton & Mathis 1989). 

We used a procedure similar to that employed for IRAS 13224-3809 and IH 0707-495 as 
described in Section 3.1 to investigate the profile of the C IV line. We first identified and subtracted 
the local continuum, and excluded portions of the profile modified by Galactic absorption lines. 
We require a smoothed profile to study the shape of the line. We used two different methods 
depending on the signal-to-noise ratio to obtain the smoothed profile. When the signal-to-noise 
ratio was relatively high, we fit the continuum-subtracted spectrum in the region of the emission 
line with a spline function as described in Section 2.3.1. When the signal to noise is lower, the 
multinode spline could not be constrained, so instead we fit the profile with several Gaussians or 
Lorenzians, and used the model as our smoothed profile. In a few objects (Ark 564, WPVS 007, 
PG 1404+226), the line profiles were altered by the presence of strong associated absorption lines. 
For these objects, the line profiles were reconstructed using multiple-component fits, but there is 
no doubt that the result is uncertain, especially with regard to the equivalent widths. Therefore, 
these objects are denoted by a different symbol in plots that include C IV information. A broad 
absorption line is observed in PG 1001+054 (Brandt, Laor & Wills 2000). However, this broad 



'■http://www.stecf.org/poa/pcrel/POA_CALFOS.html 
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absorption line appears to be sufficiently detached that it does not alter the C IV line profile. 

We then measured a number of properties from the resulting C IV profile, including the fiux 
and equivalent width. We investigated several possible ways to parameterize the shape and offset of 
the C IV line. We considered using the velocity offset (e.g., Marziani et al. 1996); however, this does 
not take into account the fact that the lines are very asymmetric. We looked into using the moments 
of the line profile (mean, variance, skew and kurtosis). We also investigated using the fraction of 
the line blueward of the rest wavelength (here taken to be 1549.5 A, the average of the doublet 
wavelengths). We found that the mean wavelength is strongly correlated with the fraction of the 
emission blueward of 1549.5 A, and finally opted to use the latter as the parameter representing 
the asymmetry and blueshift of the line. We also obtained the kurtosis as a measurement of the 
peakiness of the line: a line with positive kurtosis is pointy, and a line with negative kurtosis has 
a flattened profile. 

Other high ionization lines were measured. Wc obtained a rough estimate of the flux and 
equivalent width of N V by modeling it with several Gaussians or Lorenzians. In some objects, 
the lines were narrow enough that N V could be cleanly separated from Lya, and the principal 
uncertainty in the measurements of the fiux and equivalent width is the placement of the continuum. 
In other objects, N V overlapped Lya; we estimate the uncertainty in the fiux and equivalent width 
may be as high as 40% for these objects. As discussed in Section 2.3.3 , the O IV] and Si IV lines 
near 1400A were quite difficult to separate in IRAS 13224-3809 and IH 0707-495. Thus, for these 
other objects, we measure the fiux and equivalent width of the whole 1400A feature by simply 
integrating over it after defining a local continuum. He II was also measured. In some objects, this 
line was blended with O III]AA1660.8, 1666.1; in these cases, we used a model consisting of several 
Gaussians to isolate the contribution of He II. In several objects, a significant source of uncertainty 
was determination of the continuum when the relatively weak He II line was blueshifted and broad. 
Finally, we computed ratios of the fiuxes of the high-ionization lines with that of C IV. 

The feature near 1900A is comprised primarily of the intermediate ionization lines Al III, Si III], 
C III], and Fe III UV34. This was modeled in a way similar to that described in Section 2.3.4. As 
discussed in Vestergaard & Wilkes (2001), the Fe III UV34 multiplet (1894.8A, 1914. lA, 1926.3A; 
^5— ^P°) contributes significant uncertainty to line measurements in this region. However, in all of 
the objects that we looked at, the intermediate ionization lines appeared to be narrow, symmetric, 
and centered at their rest wavelengths. We therefore assumed that all of the lines in the 1900A 
feature have the same width, and assumed that they are centered at their rest wavelength. The 
presence of Fe III UV34 generally could be ascertained by the presence of the 1926A line, which 
was discernible from the wings of C III] lines in these narrow-line objects. The gf values obtained 
from The Atomic Line list^ indicate that in the optically-thin case, the line ratios 1926:1914:1895 
should be 1:1.4:1.8. However, radiative transfer is likely to be important for these lines. In some 
cases, the spectra were fit well assuming that the three Fe III lines have equal intensity, as would 



'http://www.pa.uky.edu/~peter/atomic/ 
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be appropriate if the gas is very optically thick. However, in several objects (e.g., Ark 564; Fig. 
15), the 1914A component was much stronger than the other two; this enhancement was discussed 
previously by Vestergaard & Wilkes (2000) for I Zw 1, although these authors do not give an 
explanation for it. It turns out that the upper-level energy for this transition is 82333. 92cm~^, 
corresponding to 1214.6A, which is just l.lA from Lya. Therefore, the 1914A component of the 
Fe III UV34 feature is plausibly selectively excited by Lya. Interestingly, this seems to be almost 
the only line emitted under broad-line region conditions that can be excited from the ground state 
by Lya; many more lines, including O IA1302 for example, may be pumped by Ly/3, and Fe II 
may be pumped by Lya from low-lying metastable levels (e.g., Sigut & Pradhan 2003). Therefore, 
we also fit with a model in which the 1895 and 1926 A components are constrained to have equal 
intensities, and the 1914 A intensity is left free. In this case, the 1914 A component ranges from 1.8 
to typically 3 times the strength of the other two. The fluxes of Al III, Si III], and C III] lines, and 
the velocity width of C III] were measured, and the ratios of Al III and Si III] to C III] were formed. 
The 1900 A feature is complex; thus, there is significant uncertainty in the relative intensities of the 
Si III] and C III] lines; Al III was separated sufficiently from the the other two that uncertainty in 
this line was dominated by the continuum placement. 

We obtained the equivalent widths and velocity widths of the Mg II line after subtracting 
the UV Fe II as described in Section 2.3.6. A by-product of the Fe II subtraction procedure is a 
measure of the slope of the continuum a„, where F{\) oc A°", in the UV between about 2200A 
and 3050A (see Moore & Leighly, in prep, for details). We also measured the monochromatic 
luminosity at restframe 2500A and 1400A. Finally, we used generally nonsimultaneous X-ray data 
to estimate aox, the point-to-point power-law energy index between 2500A and 2 keV in the rest 
frame. If available, ASCA data were used because ffux at 2 kcV could be estimated better than 
from ROSAT data in which 2 keV falls at the end of the bandpass. Uncertainties in intrinsic aox 
originate in the lack of simultaneity of the observations, and reddening or absorption arising in the 
host galaxy or AGN. 

3.1. Correlations and Principle Components Analysis 

We performed a correlation analysis of the parameters discussed above, using the Spearman 
rank correlation coefficient. In a few cases, parameters could not be measured: an example is 
RX J0134.2-4258, in which no C III], Si III] and Al III were detected. We also do not use the 
velocity width of C III] for IH 0707—495 because it was much larger than the widths of other 
similar lines in that object (sec Section 2.3.4). Since we are only missing one object (three in one 
case) among 16, we do not use a survival analysis. The probabilities of accidental correlation are 
given in the correlation matrix in Table 3. 

We performed a principle components analysis using the correlation matrix. Although the 
number of objects and parameters is small, we obtain a fairly significant clustering of correlations. 
The first three eigenvectors are given in Table 4; they represent 78% of the variance in the sample. 
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We caution that our sample is small, so these eigenvectors may not survive if faced with a larger 
sample. Nevertheless, they provide a convenient organizational framework in which to discuss the 
correlations. 

3.1.1. Eigenvector 1: the C IV and other high-ionization line properties 

Eigenvector 1 is aligned along the mutually strong correlations between C IV equivalent width 
and profile shape. A number of other parameters are strongly correlated or anticorrelated with the 
asymmetry of the C IV line. These include the high-ionization properties of the equivalent width of 
He II and the N V/C IV ratio, the intermediate- ionization line properties of the Al III/C III] ratio 
and the Si III]/C III] ratio, and the continuum property aox- A number of these are shown in Fig. 
7. Another interesting correlation is the one between the N V/C IV ratio and the Al III/C III ratio; 
we defer that discussion until Section 3.1.3. 

The C IV asymmetry parameter is strongly anticorrelated with the C IV equivalent width; that 
is, the strongly blueshifted lines have the lowest equivalent width. This was previously reported in 
Leighly 2001 for this sample, and has since then been recovered from the SDSS quasars (Richards 
et al. 2002), and a similar correlation was previously noted by Wills et al. (1993). The kurtosis is 
strongly anticorrelated with the asymmetry parameter, implying that the more strongly blueshifted 
lines have a flattened rather than peaked profile. This result suggests the hypothesis that the 
narrow, symmetrical part of C IV is an additional component that increases the equivalent width 
of the whole line when present, as well as reducing the cumulative offset from the rest wavelength. 

To test this idea, we constructed a composite line consisting of the wind profile from IRAS 13224—3809, 
and a narrow and symmetric line that has a root-secant profile with a FWHM of 2000 km s~^ and 
varying normalization. This FWHM was chosen because with it we could best match the rela- 
tionship between kurtosis and asymmetry parameter (dashed line in the second panel of Fig. 7) 
when used in combination with the IRAS 13224—3809 wind profile. We assumed that the wind 
component has a fixed equivalent width, and that the narrow symmetric component contributes 
with strength varying from zero in the case of IRAS 13224—3809. We then measured the asym- 
metry parameter for these simulated line profiles, and infer an equivalent width by scaling the flux 
with that of IRAS 13224—3809. This relation is shown by the dashed line in the top panel of Fig. 
7. Surprisingly, despite the very simple assumptions, the correspondence is quite good. We note 
that at least one other factor should influence the line equivalent width: some objects have bluer 
contimia than others. In fact, there is an anticorrelation between and the asymmetry parameter, 
such that objects with asymmetrical lines have bluer continua, and therefore would have a lower 
equivalent width for the same C IV line flux. We note that this is not a clear manifestation of the 
anticorrelation between line equivalent width and luminosity (the Baldwin effect; Baldwin 1977); 
Fig. 8 shows that the objects with the lowest equivalent widths have mid-range luminosities. We 
discuss the Baldwin effect further in Section 3.4.3 
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The asymmetry parameter is anticorrelated with both continuum parameters aox and a„. 
Thus, objects with blueward asymmetric lines have relatively bluer UV continua and relatively 
lower X-ray flux. An exception is I Zw 1, which appears to have a reddened UV continuum (Smith 
et al. 1997). This dependence can be directly interpreted in terms of resonance line-driven winds in 
AGN. In such winds, the acceleration comes from resonance scattering of the UV continuum from 
the central engine or accretion disk by ions in the wind. An intrinsically blue continuum provides 
more photons to be scattered; thus a blue continuum should be associated with an asymmetric 
line profile, as is seen. Both Murray et al. (1995) and Proga, Stone & Kallman (2000) discuss the 
fact that strong soft X-ray emission tends to overionize the ions required for resonance scattering; 
without these ions, the wind is suppressed. Thus, lower X-ray flux should be associated with the 
presence of the asymmetrical line profile, as is observed. 

We note that very few of the UV and X-ray observations were coordinated, so relative variability 
between UV and X-ray observations adds noise to aox', however, such noise would tend to decorrelate 
the parameters, rather than cause a false correlation. Also, a few objects may be reddened (e.g.. 
Ark 564; Crenshaw et al. 2002); this would artificially increase aox- Also, PG 1001-1-054 has a 
broad absorption line (Brandt, Laor & Wills 2000), and a very steep ctorr, which may imply that 
the X-rays are absorbed. Nevertheless, when the suspect objects are removed, the trend remains 
{P = 1.9%), tentatively supporting our hypothesis. Furthermore, we note that our hypothesis is 
supported by the emission-line properties of two objects for which we do have simultaneous UV and 
X-ray observations (RE 1034-39: Casebeer & Leighly 2004; PHL 1811: Leighly et al. in prep.). 

Also associated with the asymmetry parameter and C IV equivalent width is the equivalent 
width of He II; objects with blueshifted C IV have weak He II lines. This result has two interpre- 
tations. It may be a consequence of the fact that in any particular object, if one line has a high 
equivalent width, many of the other lines also have a high equivalent widths, possibly a result of 
a deficient continuum or large covering fraction, for example. However, there is another interpre- 
tation for the anticorrelation between the asymmetry parameter and the He II equivalent width. 
A weak He II line is generally considered to be smoking-gun evidence for a continuum weak in 
soft X-rays; if this is the factor driving the correlation, it would also support the presence of a 
resonance-line driven wind, as mentioned above. 

Interestingly, there is no similar correlation with the equivalent width of the stronger high- 
ionization line, N V. There is an weak correlation between N V equivalent width and a^, such that 
objects with lower N V equivalent widths have bhicr spectra. This could be because bluer continua 
yield lower equivalent widths for the same line flux. Stronger correlations are seen with the ratio 
of N V to C IV. Some of these may be related to C IV correlations, but others are more strongly 
correlated with the ratio than with the C IV alone. The anticorrelation between He II and the ratio 
of N V to C IV may support the idea that objects having a wind may have N V selectively excited 
by Lya; this issue is discussed further in Paper II. Several correlations are found between this ratio 
and properties of the intermediate-ionization lines; these will be discussed in Section 3.1.3. 
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3.1.2. Luminosity/Equivalent Width Correlations 

Wc found a strong anticorrelation between the C IV equivalent width and the profile asym- 
metry. It has been known for years that there exists generally an anticorrelation between AGN 
emission line equivalent widths and luminosity; this is known as the Balwin effect (Baldwin 1977). 
The physical origin of the Baldwin effect has remained a puzzle for many years; however, recent 
work suggests that it is a result of UV/X-ray continuum becoming softer as the luminosity increases 
(Dietrich ct al. 2002). In this sample, the correlations reported in Table 3 between line equivalent 
width and the luminosity at 2500A are weak; it may be that the correlations are not stronger 
because we probe only a limited range of luminosity. 

In Fig. 8 we compare the relationship between XLx and equivalent width for our sample with 
the regression obtained for a much larger sample by Dietrich et al. (2002). We note that for the 
purpose of this plot, the luminosities were calculated assuming the cosmology used by Dietrich et 
al. (2002): Hq = 65kms~^Mpc, = 0.3, and U\ = 0.7. The most remarkable feature of this 
plot is that, although there is still an anticorrelation between equivalent width and luminosity for 
NLSls, the relationship for many of the emission lines is offset from that of the average quasar; for 
a given luminosity, NLSls have lower equivalent-width emission lines. This peculiar property of 
NLSls was previously reported by Wilkes et al. (1999). Interestingly, N V is the exception to the 
trend; the N V equivalent widths from NLSls are consistent with those from the average quasar. 

IH 0707-495 and IRAS 13224-3809 have the smallest C IV equivalent widths in the sample; 
however, they have median luminosity, so their small equivalent width is not a consequence of the 
Baldwin effect. Interestingly, for some of the lines, the equivalent widths from these two objects is 
consistent with those of the other NLSls; for other lines, they are markedly lower. 

3.1.3. Intermediate-ionization Line Properties 

Eigenvector 2 is aligned along the equivalent widths of the Si III] and Al III; the equivalent 
widths of these lines arc strongly correlated. Mutual correlations are found between the ratio of 
C III] to C IV, the equivalent widths of Al III and Si III], and the ratios of Si III] and Al III to C III]. 
These interdependencies are illustrated in Fig. 9. 

What is the origin of these interdependencies? We discuss this qualitatively here, and provide 
a more quantitative discussion, including photoionization modeling, in Paper II. First, it is possible 
that the trends reflect a difference in density of the line-emitting region. Because of its low critical 
density of 3.4 x 10*^ cm ^, C III] is expected to decrease with respect to Si III] (critical density of 
1.04 X 10^^ cm~^) and Al III (a permitted line) as the density increases. Then, as density increases, 
Si III] and Al III are expected to take over the cooling, resulting in an increase in both their 
equivalent widths and their ratios with respect to C III]. However, an increase in density would 
predict a decrease in the C in]/C IV ratio, since the critical density of C IV is 2.06 x 10^^ cm^^. This 
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would predict an anticorrelation between C III] /C IV and Si III] / C III] , rather than the correlation 
that we observe. 

Another possible origin of these trends is a variation in the ionization parameter. A decrease 
in ionization parameter predicts an increase in the ratio of C III] to C IV. It may conceivably 
produce larger equivalent widths of Si III] and Al III, because their lower ionization potential allows 
them to be produced under lower ionization conditions, leading to larger ratios of Si III]/C III] and 
Al III/C III]. 

There is another possible explanation that we discuss in much greater detail in Paper II; 
photoionization modeling is presented there that supports this explanation. If we assume that the 
objects that have strongly blueshifted C IV lines are characterized by the presence of a wind, it 
is possible that the continuum passes through the wind, or is "filtered" by the wind, before it 
illuminates the intermediate-ionization line-emitting region^. After passing through the wind, the 
continuum will be somewhat weaker in the hydrogen continuum and will have very few photons in 
the helium continuum. Such a continuum is effectively softer overall, and therefore lower-ionization 
species should dominate the cooling. As the continuum becomes softer, we see a decrease in the 
ratio of C III] to C IV, because the continuum is no longer hard enough to ionize 0+2 efficiently. 
Al III and Si III] then increase relative to C III] because Al-l-2 and Si+2 have lower ionization 
potentials than C+2. 

Some of the intermediate-ionization emission line properties are correlated with properties of 
the high- ionization lines. There are strong correlations between Al III/C III] and Si III/C III] ratios 
and the asymmetry parameter (Fig. 7). This perhaps also supports the concept of filtering outlined 
above, since that can be interpreted as a shift of the cooling in the intermediate-line emitting region 
to ions with lower ionization potentials when a wind is present. We also see a correlation between 
the intermediate-ionization line ratios and the N V/C IV ratio that is particularly strong for the 
Al III ratio. Because N V overlaps Lya in the windy objects, it is possible that N V is selectively 
excited by Lya in these objects. Thus, this correlation could be interpreted as additional support 
for the filtering scenario outlined above. 

3.1.4. The Line Velocity Widths 

The velocity widths appear in the third eigenvector. The correlations between H/3, Mg II and 
C III] velocity widths are not as strong as might be expected, as seen in Fig. 10. The scatter 
possibly originates in nonuniform measurement of the H/3 widths; many of these values were taken 



®We differentiate between a "shielded" continuum, which is assumed to have been transmitted through highly 

ionized gas (e.g., Murray ct al. 1995), and a "filtered" continuum, which is assumed to have been transmitted through 
the wind while ionizing and exciting it before illuminating the disk and producing the observed intermediate- and 
low-ionization lines. This point is expanded upon in Paper II. 
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from the literature. Also, the width of C III] is difficult to measure, as it is blended with Si III] and 
Fe III UV 34. We do observe a correlation between the velocity width of Mg II and the luminosity 
at 2500A. This correlation is expected, as objects with higher luminosities should have larger black 
holes and correspondingly BLR velocities (Laor 1998); we discuss this correlation in a larger sample 
in Leighly Sz Moore 2004 (in prep.). We also see anticorrelations between the C III]/C IV ratio 
and H/3 velocity width, and the He II/C IV ratio and Mg II velocity width. We have no ready 
explanation for these correlations. 



4. Summary 

• We present a detailed analysis of the UV spectra from two Narrow-line Seyfert 1 galaxies, 
IRAS 13224-3809 and IH 0707-495. We find that their continua are as blue as that of the av- 
erage quasar. We observe that the high-ionization emission lines (including N V and C IV) are 
broad (FWHM?a 5000 kms~^), and strongly blueshifted, peaking at around 2500 kms"^ and 
extending up to almost ^ 10, 000 km s~^. In contrast, the intermediate- and low-ionization 
lines (e.g., C III] and Mg II) arc narrow (FWHM 1000 1900 km s"^) and centered at the 
rest wavelength. Si III] is prominent, and other low ionization lines (e.g., Fc II and Si II) 
are strong. Based on these observations, the working model that wc adopt considers that the 
blueshifted high-ionization lines come from a wind that is moving toward us, with the receding 
side blocked by the optically thick accretion disk, and the intermediate- and low-ionization 
lines are emitted in the accretion disk atmosphere or low- velocity base of the wind. 

• The strongly blueshifted C IV profile suggests that it is dominated by emission in the wind. 
We used the C IV profile to develop a template for the wind. We then used this template, 
plus a narrow and symmetric component representing the disk emission to model the other 
bright emission lines. We inferred that the high-ionization lines N V and He II are also 
dominated by wind emission, and a part of Lya is emitted in the wind. The intcrmediate- 
and low-ionization lines Al III, Si III], C III], and Mg II are dominated by disk emission, and 
a part of Lya is disk emission as well. The 1400A feature, comprised of O IV] and Si IV was 
difficult to model; however, it appears to include both disk and wind emission. 

• IRAS 13224-3809 and IH 0707-495 have distinctive X-ray properties among NLSls; in 
order to determine whether their distinctive properties carry over to the UV, we analyzed 
HST archival spectra from 14 other NLSls with a range of two orders of magnitude in UV 
luminosity. We find indeed that these two objects are extreme in this sample in the following 
properties: strongly blueshifted C IV line, the low equivalent widths of many of the lines, 
particularly C IV and He II, high C III]/C IV, Si III]/C III], Al III/C III], and N V/C IV 
ratios, steep aox, and blue UV continuum. Correlation analysis finds a number of strong 
correlations. An anticorrelation between C IV asymmetry and equivalent width suggests 
that the line is generally composed of a highly asymmetric wind component and a narrower 
symmetric component. The anticorrelation between C IV asymmetry and aox and and 
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with He II suggests that UV-strong and X-ray weak continua may be associated with a wind, 
as would be expected if the acceleration mechanism is radiative-line driving. The dominance 
of Al III and Si III] over C III] suggests possibly that the continuum is transmitted through 
the wind before it illuminates the intermediate-ionization line emitting gas. 

• Paper II investigates the physical conditions in the line-emitting gas through Cloudy modeling 
of the disk and wind emission lines. The results are used, in combination with a toy dynamical 
model for the wind, to estimate the distance of the wind from the continuum emitting region. 
Further discussion of the results is given, including a comparison with previous results. 
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1000 2000 5000 

wavelength (A) 

Fig. 1.— The broadband continua of IH 0707-495 (top) and IRAS 13224-3809 (middle) compared 
with that of a composite quasar spectrum (bottom; Francis et al. 1991). The dashed hne is a power 

law with slope of —7/3, which is the slope predicted by the standard thin accretion disk model. 
The dotted line is a power law with slope —0.32, which is the median in the LBQS (Francis et al. 
1991). 
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Fig. 2.— The HST STIS spectra from IRAS 13224-3809 and IH 0707-495, shifted to the rest 
frame and dereddened as described in the text. For clarity, spectra have been renormahzed to have 
the same flux, and offset from one another. The rest frame wavelengths of strong emission lines 
frequently seen in AGNs are labeled, excluding the many Fe II multiplets. All identified absorption 
lines are labeled; they all originate in our Galaxy. 
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Fig. 3. — The rescaled representative high-ionization line C IV superimposed on the representative 
low- ionization line Mg II as a function of velocity, for our two NLSls and the average quasar (Francis 
et al. 1991). The average quasar high-ionization line is slightly broader and slightly blucshiftcd 
compared with the low-ionization line. In contrast, the low-ionization lines in NLSls are much 
narrower and the high-ionization lines are strongly blueshifted. 
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Fig. 4. — Line models for the C IV region, the Lycc/N V region, and the Si IV/0 IV] region. The 

dark grey dashed Hnes show the components and the hght grey vertical lines show the positions 
of the lines in the rest frame. Resonance lines were assumed to be optically thick and thus their 
ratios were assumed to be 1:1. O IV] was assumed to be optically thin (see text). 
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Fig. 5. — Line models for the Al Ill/Si III]/C III] region, the Mg II region, and the H/3/[0 III] 
region. The dark grey dashed hnes show the components and the hght grey vertical lines show the 
positions of the lines in the rest frame. 
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Fig. 6. — Model for the 1900A feature in Ark 564. Diamonds mark the rest wavelength positions of 
the emission lines. The very narrow lines allow the pumped Fe IIIA1914 feature to be unambiguously 
detected. 
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Fig. 7. — Various parameters as a function of the asymmetry parameter. Stars mark 
IRAS 13224—3809 and IH 0707—495, and open circles mark objects with significant C IV associated 
absorption lines; their profiles and equivalent widths are uncertain because of assumptions made to 
reconstruct the profiles (see Section 3). Equivalent widths (EW) arc given in units of Angstroms. 
aox {F{E) oc is defined between 2500A and 2 keV in the rest frame. a„ {F{X) oc A"") is 

obtained as a byproduct of the UV Fe II subtraction, and is defined between approximately 2200 
A and 3050 A. Top and second panels: the dashed line shows the predicted relationship between 
equivalent width and kurtosis and asymmetry parameter for a two-component line comprised of a 
wind and symmetric component (see text). 
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Fig. 8. — Measured equivalent widths as a function of XLx at 1450A. Symbols have the same 
meaning as in Fig. 7. Dashed lines are the Baldwin effect regressions computed by Dietrich et al. 
2002. Note that for the purpose of this plot, luminosities were computed using the cosmological 
parameters assumed by Dietrich et al. (2002). 
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Fig. 9. — Intermediate-ionization line equivalent widths as a function of the ratio of C III] to C IV 
and the ratio of Si III] to C III]. Symbols have the same meaning as in Fig. 7. 
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Fig. 10. — Velocity widths of C II] and Mg II as a function of H/3 velocity width. Stars denote 
IRAS 13224-3809 and IH 0707-495. 
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Table 1. Observing log 



Target 


Date 


Spectrometer 


Observed Wavelength 


Resolution 


Aperture^ 


Exposure 








(A) 


(A) 


(arc seconds) 


(seconds) 


IRAS 13224-3809 


1999-06-05 


HST STIS G230L 


1568-3184 


3.2 


0.5 X 0.24 


4457 




1999-06-05 


HST STIS G140L 


1140-1730 


1.2 


0.5 X 0.24 


8188 




1999-06-20 


CTIO 4m 


3870-7570 


3.0 


1x4 


3600 


IH 0707-495 


1999-02-16 


HST STIS G230L 


1568-3184 


3.2 


0.5 X 0.24 


2196 




1999-02-16 


HST STIS G140L 


1140-1730 


1.2 


0.5 X 0.24 


5376 




1998-01-04 


CTIO 1.5m 


3500-6940 


5.5 


1.8 X 7 


1800 



Slit width X extraction aperture in the spatial direction 



Table 2. Emission Line Measurements 



Emission line^ 




IRAS 13224-3809 






IH 0707-495 






Average 


EWs (A) 




Flux 


Equivalent Width 


Width 


Flux 


Equivalent Width 


Width 


Francis 


Zheng'' Brotherton 




(xlO~l*erg cm~2 


s-i) (A) 


(kms-i) 


(xlO~^*erg cm~2 J 


i-i) (A) 


(kms-i) 








Lya+N V feature 


59.4 


42 




166 


45.9 




52 




87 


Lya (broad) 


19.9 


14.0 




65.3 


17.7 






85 




Lya (narrow) 


8.3 


6.0 


1160 


20.2 


5.6 


1135 








N V (broad) 


27.9 


20.5 




57.8 


16.4 






10 




1400 A feature 


11.4 


9.6 




28.3 


9.9 




10 


7.3 


8 


Si IV (broad) 


5.5 < 7.2« 


4.6 < 6.1'' 




10.1 < 15.9'= 


3.5 < 5.6'= 










Si IV (narrow) 


0.75 < 3.7"^ 


0.63 < 3.2'= 




4.5 < 10.6"= 


1.6 < 3.7'= 










O IV (broad) 


3.7 < 10.2"^ 


3.1 < 8.6'= 




9.4 < 22.4<= 


3.3 < 7.9'= 










O iV (narrow) 


1.3 < 3.4° 


1.1 < 2.9'= 




2.3 < 7.6° 


0.8 < 2.7° 










C IV (broad) 


14.6 


12.8 




31.3 


11.6 




37 


59 


33 


He II (broad) 


3.1 


2.8 




8.3 


3.2 




12<* 


3.9 


7.0=^ 


1900 A feature 


17.0 


19.0 




42.5 


18.3 




22 




1^ 1 


Al III 


3.5 


3.7 


1835'' 


8.6 


3.7 


3500" 




3.5 


CO 


Si III] 


6.5 


7.2 


1835*= 


15.7 


6.8 


3500== 




3.5 


00 


ClII] 


4.7 


5.3 


1835'= 


11.9 


5.2 


3500-= 




17.0 


1 


Mg II 


6.7 


17.3 


866'' 


11.6 


10.7 


1050* 


50 


64 


34 


C IIA1335 


3.5 


2.9 


1700 


8.2 


2.7 


2100 








C IIA2327 


1.5 


2.7 
















N IIIA1750 


1.0 


1.2 




4.4 


2.0 











""Broad" refers to lines fit with the C IV template presumably coming from the wind, while "narrow" refers to lines symmetric at zero velocity presumably coming from the 
disk or low velocity base of the wind. 

''Radio-quiet quasars. 

^Best value and upper limit; see text for details. 
'^He IIA1640+O III]A1663 
''Widths constrained to be equal in fitting. 
'^Width of one component of the doublet. 



Table 3. Correlation Matrix 



Property 




C IV Properties 




High-ionization Line 


Properties 




Intermediate-ionization Line Properties 


C IV EW 


Asymmetry 


Kurtosis N V EW 


N V/O 1 V 


TT TT 77(11 r 

He 11 


He II/C IV 


1400A EW 


C III] EW C III]/C IV C III] FWHM 


cIVew 




-5.2 X 10-3 


0.10 


—0.40 


1.4 x 10-3 






2.1 ••• -4.7 


Asymmetry 


-5.2 X 10-3 




-2.8x10-5 


0.16 


-0.023 






-1.7 ■•• 5.2 


Kurtosis 


0.10 


-2.8 X 10-5 




-0.24 


0.057 






2.8 • ■ ■ -3.4 


N V EW 


















N V/C IV 


-0.40 


0.16 


-0.24 




—0.51 






-0.069 2.2 


He II EW 


1.4 X 10-3 


-0.023 


0.057 


-0.51 








0.16 


He II/C IV 


















1400A EW 










3.9 






1.8 


C III] EW 


2.1 


-1.7 


2.8 


-0.069 


0.16 








C III]/C IV 








2.2 




2.0 






C iii] FWHM 


—4.7 


5.2 


—3.4 












Al III EW 
















1 


Al III/C III] 


-1.4 


0.078 


-0.022 


8.0 x 10-3 


-0.26 






— 1.3 ••• 


Si III] EW 


















Si III]/C III] 




0.86 


-0.22 


0.30 








17 1 ... 


Mg II EW 


















Mg II FWHM 












-0.45 






H/3 FWHM 












-4.2 




-0.29 


Otox 


0.028 


-0.50 


0.48 




0.043 


4.7 




-0.87 


Oiu 


1.9 


-2.7 


5.5 2.5 


-4.4 


0.45 




0.30 


3.1 


^2500 
















-4.4 



Note. — 



Percentage probability that the observed correlation is accidental are listed when the value is less than 6 (i.e., greater than 94% confidence that correlation is real). 



Table 3. Correlation Matrix (cont.) 



Intermediate-ionization Line Properties (cont.) Low-ionization Line Properties Continuum Properties 



Property 


Al III EW 


Al III/C III] 


Si III] EW 


Si III]/C III] 


Mg II EW 


Mg II FWHM 


H/3 FWHM 


OLox 




^-2500 


C IV EW 




-1.4 












0.028 


1.9 




Asymmetry 




0.078 




0.86 








—0.50 


—2.7 




Kurtosis 




—0.022 




—0.22 








0.48 


5.5 




N V E W 


















2.5 




N V/C 1 V 




8.0 X 10 




0.30 










—4.4 




ne 11 r!j vv 




— xy.zXi 












U.U^to 


U.40 




TT TT //^ TA J 

He ii/C i V 












—0.45 


—4.2 


4.7 






1400A EW 






2.3 




2.5 








0.30 




C III] EW 




-1.3 






1.1 








3.1 




C III]/C IV 








1.7 






-0.29 






-4.4 


C III] FWHM 
















-0.87 






Al III EW 






1.3 X lO"'' 


0.79 


2.4 












Al III/C III] 








0.010 








-4.3 






Si III] EW 


1.3 X 10"" 






2.4 


0.28 












Si III]/C III] 


0.79 


0.010 


2.4 
















Mg II EW 


2.4 




0.28 
















Mg II FWHM 




















0.40 


H^ FWHM 




















5.9 


O.OX 




-4.3 


















Oil 


































0. 10 


5.9 









Note. — Percentage probability that the observed correlation is accidental are listed when the value is less than 6 (i.e., greater than 94% confidence that 
correlation is real). 
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Table 4. Principal Components Analysis Results 



Parameter 


Eigenvector 1 


Eigenvector 2 


Eigenvector 3 


CIV EW 


0.32 


0.02 


0.01 


Asymmetry 


-0.32 


0.06 


0.01 


Kurtosis 


0.32 


-0.08 


-0.09 


N V EW 


0.02 


0.21 


0.27 


N V/C IV 


-0.29 


0.05 


-0.24 


He II EW 


0.33 


0.07 


-0.05 


He II/C IV 


0.06 


0.15 


-0.42 


1400A EW 


0.17 


0.26 


0.18 


C III] EW 


0.26 


0.16 


0.14 


C III]/C IV 


-0.11 


0.27 


-0.33 


C III] FWHM 


-0.21 


0.06 


0.31 


Al III EW 


-0.06 


0.40 


0.16 


Al III/C III] 


-0.30 


0.16 


-0.02 


Si III] EW 


-0.01 


0.42 


0.15 


Si III]/C III] 


-0.23 


0.32 


-0.02 


Mg II EW 


0.12 


0.26 


0.27 


Mg II] FWHM 


-0.09 


-0.23 


0.30 


H^j FWHM 


0.06 


-0.23 


0.34 




0.28 


-0.02 


-0.20 




0.24 


0.19 


0.12 


L25OO 


-0.17 


-0.27 


0.20 


Cumulative Fraction of Variance 


0.38 


0.60 


0.78 



